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Abstract 


In a Sikivie-type cosmic-axion detector, both the width and position of the microwave 
signal due to axion-photon conversion depend upon the motions of the earth. Due to 
the orbital and rotational motions of the earth they will be modulated with periods of 
1 sidereal day and 1 sidereal year, with amplitudes of about 0.1% and 5% respectively. 
Because of the intrinsically-high energy resolution of Sikivie-type detectors such periodic 
variations should be detecta ble. Such modulations wou ld not only aid in confirming the 
detection of cosmic axions, but, if found, would also provide important information about 
the distribution of axions in the halo. 
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The axion is a very attractive dark-matter candidate. Moreover, at present here are 
several ongoing experiments 1 to detect cosmic axions in our halo by the process proposed 
by Sikivie: 2 Resonant conversion of axions into microwave photons in the presence of a 
strong magnetic field. In addition, a second generation experiment that is about 300 times 
more sensitive is being planned. 3 The characteristic energy of the photons produced in a 
Sikivie detector is = m a c 2 + m a i > 2 /2 where m a is the axion mass and v is the velocity 
dispersion of axions in the halo (we will be much more precise later). For closure density, 
m 0 is expected to be in the range 1(T 6 eV to 10" 4 eV (Refs. 4); the velocity dispersion of 
our halo is about 300 km sec" 1 ~ 10" 3 c, so that m a v 2 / 2 ~ l(r 6 m a . The frequency of the 
axion microwave line is 2.4 GHz(m a /10 5 eV), while the width corresponds to a frequency 
of about (m a /10 5 eV)kHz. Because resonant- cavity experiments have very good energy 
resolution Q s in excess of 10 5 both the position and the shape of the microwave line 
are potentially observable. 

Since the earth is in motion with respect to the halo — its motion as part of the solar 
system, its orbital motion, and its rotational motion — the position and width of the mi- 
crowave line expected from axion-photon conversion will depend upon these velocities and 
will be periodically modulated. (We will neglect the very small motion of the earth about 
the earth-moon barycenter, v ~ 1.3 x 10- 2 km sec -1 and period 27.32 sidereal days.) The 
periodic modulations due to the earth’s orbital and rotational motions may prove crucial in 
confirming the detection of cosmic axions and, if detected, in determining the distribution 
of axions in the halo. 

It is usually assumed that the distribution of matter in the halo can be described 
by an isothermal phase-space distribution in the rest frame of the galaxy— not so much 
because of any hard evidence, rather for simplicity and the lack of a more compelling 
model. In fact, there is essentially no evidence to preclude the possibility that the halo is 
rotating with respect to the galaxy. In any case, such a phase-space distribution is a very 
reasonable starting point and will serve to illustrate the periodic effects of interest. In the 
non-rotating, inertial rest frame of the galaxy (in astrophysical parlance, the Fundamental 
Standard of Rest) the local phase- space distribution of axions is given by 


fd 3 v = n 0 



3/2 

exp(— 0v 2 )d 3 v, 


( 1 ) 


where n 0 -3x 10 13 cm 3 /(m a /10 5 eV) is the local number density of axions (to within a 
factor of about 2; see Turner in Ref. 4) and h = k B = c = 1 throughout. In all that follows 
we will ignore the factor of no, so that J fd 3 v — 1; we will also be cavalier in not relabeling 
the distribution function / when we change its argument. The quantity /? corresponds to 


m a /2T where T is the “halo temperature” and determines the velocity dispersion of the 
halo. It is very simple to carry out the angular integrations and write the distribution in 
its familiar Maxwellian form: 

fdv = 47 r v 2 exp (~/3v 2 )dv. (2) 

From this expression one can compute (u 2 ) and (v 4 ): 

< v 2 > = J v 2 fd>v = \r' <» 4 > = / = jr 2 -, 

based upon the dynamics of our galaxy 5 v = (v 2 ) 1 / 2 has been determined to be about 
270 km sec -1 , which implies that = 1.5u -2 ~ (220 km sec -1 ) -2 . The uncertainty in the 
quantity 0 is probably 10% or so. 

(Because the galaxy has an escape velocity, thought to be about 650 km sec -1 , the 
axion distribution should be truncated at this velocity. For our purposes this fact is not 
significant. In addition, Griest 6 has pointed out that owing to the gravitational effects of 
the sun there is a minimum axion velocity, the escape velocity from the solar system, locally 
about 42 km sec -1 , and a small additional modulation— due to focusing by the sun — that 
is orthogonal to the one discussed here. We will also ignore the small gravitational effects 
of the sun.) 

The earth is not at rest with respect to the galaxy; its velocity can be written as 

V£ = V S + Vo + \R , 

where Vs is the velocity of the sun with respect to the galaxy, Vo is the orbital velocity of 
the earth around the sun, and Vr is the rotational velocity of the position on earth where 
the axion detector is. To orient one, the velocity of the sun has a magnitude of about 
230 km sec -1 and is in a direction that is about o = 60° above (north of) the ecliptic 
(plane of the earth’s orbit); when projected down onto the ecliptic, it intersects the earth’s 
orbit at a point that is about 19 days before the autumnal equinox, i.e., near the earth’s 
position on 4 September. (The uncertainties in these numbers correspond to a few days, a 
few degrees and about 20 km sec -1 ; see Refs. 5) The angle between Vs and the equatorial 
plane is about x /> = 47°, and at local midnight, 19 days before the summer solstice (2 June) 
the projection of vs onto the equatorial plane and \r are parallel. The magnitude of the 
earth’s mean orbital velocity is t >o — 29.8 km sec -1 , while that of its rotational velocity at 
the equator is vr = 0.465 km sec -1 . Of course at latitude /, the rotational speed is smaller 
by a factor of cos 1. 
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What is of interest to the experimenter on earth are the velocities of halo axions with 
respect to the experiment. Let v be the velocity of an axion in the galactic rest frame, 
then its velocity with respect to the experiment is 


V a = V - V E . 

In the earth frame, the phase-space distribution of axions is given by 

f d * v ° = (“) e M~Pv 2 a )d 3 v a , (3a) 

/ n\ 1/2 

fdv a = 2 J exp(-/3v 2 - 0v E ) smh(20v E v a )dv a , (36) 

where of course v 2 a = v 2 - 2v E ■ v + v|, and Eq. (3b) foUows from Eq. (3a) by carrying 
out the angular integrations. The square of the earth’s velocity is 

v| = (v s + v 0 - l-v*) 2 , (4a) 

v £ - v s + 2 ~^o ■ v s + 2 ~VR ' , ( 46 ) 

where in Eq. (4b) we have dropped all terms that are quadratic in v 0 and v R since 
v 0 ^ 0.13i>s and vr < 0.002us. (Formally, the v 2 0 term is larger than the term involving 
v s v R ] however, we are interested only in the periodic terms, and v| is constant. If one 
wishes to, it is straightforward to retain all the terms in v^..) 

It is of some use to rewrite the distribution function in terms of the axion kinetic 
energy, E = m a v 2 a /2: 


fdE = 2 


1/2 


dE 

m a v E 


exp(-/?4 - 2/3E/m a )sinh[2^(2E/m a ) 1 / 2 i; E J; (5a) 


fdu = 



1/2 


du 


— exp[— 1.5(r 2 + u)] sinh(3r v / u); 


(56) 


where r — vs/v ~ 0.85 and the kinetic energy of the axion E = ( m a v 2 /2)u , which 
corresponds to a frequency of 980 Hz(m a /10- 5 eV)ti. The distribution function for v E = 0, 
i.e., that in the rest frame of the galaxy, as well as those for 2 June ( v E is maximum) and 
for 2 December (u^ js minimum) are shown in Fig. 1. The large effect of the motion of the 
solar system on the distribution function is apparent, as is the smaller effect of the earth’s 
orbital motion (no attempt was made to show the even smaller effect of the rotational 
motion of the earth). 
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From the distribution function we may compu te" any and all quantities of interest. 
Moreover, we can use the distribution function in any of its forms, cf. Eqs. (1), (2), (3a), 
(3b), and (5). It will be most convenient to use th e distribu tion function in the rest frame 
of the galaxy, Eq. (2). Remembering that v\ = v 2 — 2v • ve + v E it is simple to compute 
(u 2 ) and (vl) (and higher moments if wanted): 


( v l ) ~ j < 


fd\ = v\ + -r l = 4 + * 2 ; 


(vt) =Jv 


15 


fd v = —3 1 1 + -&v E + — 0 v E 

1 + ^ V E 


15-4 


. V 


, n . E 9 v 4 

1 + 2rj + 

v 1 15 v* 


E 


((»? - K’» 2 > = \r 2 


_ 2 _, 

3 U 


1 + 2^ 


(ба) 

( бб ) 

(6c) 


Substituting' in numbers, and for the moment neglecting the small contributions of the 
orbital and rotational velocities to v E , we find that 


v a = (vl) 1 / 2 ~ 355kmsec \ 


‘ 2 , a = (( u a - ( v l)) 2 ) l/2 - (305 km sec v ) 2 , 


~ 7 x 10- 7 , 


m ot’L/ 2 


m a c‘ 


maC 4 


1.9 x 10 6 


Note that because of the motion of the solar system relative to the galaxy the tths velocity 
of axions as seen on earth is about 355 km sec -1 , compared to about 270 km sec -1 as 
seen in the galactic rest frame. Since the energy of the axion-microwave line is about m a c 2 
and its width is rn a v 2 dls / 2, the final expression corresponds to 1 /Q a — the axion line is very 
narrow indeed. 

Now consider the average axion kinetic energy ( E a = m a u 2 / 2) and the dispersion in 
the axion kinetic energy (A E a = m a v dts / 2), taking into account the orbital and rotational 

motions of the earth: 

2r 2 

— ( — cos a cos uot + 

vs 


&-=£( 1+r 1 ) 


1 + 


1 + r 2 


A E a = 


v O . , V R 7 / , 

cos a cos uot H cos l cos xp cos uri 


VS 


m a v 


(7a) 


2(1 + 2r 2 )l 1/2 


1 + 


4r 2 


3 

v o . , v R 

— cos a cos uot H cos l cos xp cos uri 

Vs Vs 


]) 


1/2 


(76) 


1 + 2r 2 

where = 27r/sidereal year = 1.991 xlO -7 rad sec -1 is the earth’s orbital angular velocity, 
ur = 2 7r / sidereal day = 7.292 x 10 -5 rad sec -1 is the earth’s rotational angular velocity, 
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and time t is measured from local midnight on 2 June. The prefactors in expressions (7a, b), 
expressed as frequencies, are: m a u 2 (l + r 2 )/ 2 ~ 1.70kHz(m o /10 _5 eV) and m a v 2 [ 2(1 + 
2r 2 )/3] 1 / 2 /2 ~ 1.25kHz(m a /10- 5 eV). 

The (relative) amplitude of the sinusoidal variation of E a is: 2r 2 {yo/vs} cosct/(l -(- 
r 2 ) ~ 5.5% (orbital) and 2r 2 (v R /v s ) cos / cos 0/(1 + r 2 ) ~ 0.093% (rotational for a lat- 
itude of 37.5°). The periods are 1 sidereal year (3.1558 x 10 7 sec) and 1 sidereal day 
(86,164.091 sec) respectively. The (relative) amplitude of the sinusoidal variation of A E a is: 
2r 2 {v 0 /v s )cosa/{l + 2r 2 ) ~ 3.9% (orbital) and 2r 2 (v R /v s ) cosip cosl/(l + 2r 2 ) ~ 0.066% 
(rotational for a latitude of 37.5°). The maxima of the orbital variations occur about 19 
days before the summer solstice (2 June) and the minimum about 19 days before the winter 
solstice (2 December). 

(Because the earth’s orbit is not circular, orbital eccentricity e = 0.016750, there are 
higher harmonics in the moduations of E a and AE a , proportional to e N_1 cos Nu> 0 t- To 
order e, the orbital velocity of the earth is vq = au>o[~(sinu>c>t + e sin 2 u>ot)x + (cos u>ot + 
e cos2u>o<)y], where a = 1AU = 1.496 x 10 13 cm is the semi-major axis of the earth’s 
orbit, time is measured from perihelion (approximately 3 January), and at perihelion the 
earth is located on the positive x-axis. 7 The sinusoidal variation of E a at frequency 2u>q is 
2er iyo/vs) cos a cos 2u>ot' /{\ + r 2 ) and that of A E a is 2er 2 (t>o/i>s) cos o cos 2u>ot' / {]- + 
2r ), each smaller than that at frequency ujq by a factor of e ~ 0.017. Here time t' is 
measured since 20 March.) 

Note that the amplitudes of the periodic variations of E a and A E a depend upon the 
astronomical parameters v 0 /v s , v R /v s , a, ip, l, and r = vs/v. Moreover, the functional 
dependences of the modulation factors for E a and A E a are different. Thus, if these effects 
are measured, they will prove extremely useful in probing the distribution of axions in 
the halo. For example, the periodic orbital variations of E a and A E a differ by a factor 
of (1 + 2r 2 )/(l + r 2 ), a fact which could be used to infer v and therefore /?. In addition, 
the phases of the modulations can be used to infer the direction of the sun’s velocity, as 
projected onto the ecliptic. 

The importance of the yearly modulation due to the earth’s orbital motion in bolo- 
metric dark matter detection experiments has been previously emphasized. 8 There, the 
modulation is typically expected to be of order 10% and is crucial to identifying the 
signal in that case energy-deposition rate — as being due to dark matter interactions in 
the detector. Whether such detectors can achieve the required sensitivity and energy reso- 
lution to observe periodic variations remains to be seen. While the yearly modulations in 
an axion detector are slightly smaller, owing to the excellent intrinsic energy resolution of 


6 



V 


Sikivie-type detectors these variations should be easily detectable, and can be used both 
to confirm that the signal is associated with cosmic axions and to probe the halo distribu- 
tion of axions. 9 It may even be possible to detect the daily variations due to the earth’s 
rotation. 

This work was supported in part by the DOE through grant FG02 90ER40560 at 
Chicago. 


References 

1. S. DePanfilis et al., Phys. Rev. Lett. 59, 839 (1987); Phys. Rev. D 40, 3153 (1989); 
C. Hagmann, P. Sikivie, N.S. Sullivan, and D.B. Tanner, Phys. Rev. Lett., in press 
(1990); Rev. Sci. lustrum., in press (1990). Also see, Proceedings of the Workshop on 
Cosmic Axions, eds. C. Jones and A. Melissinos (WSPC, Singapore, 1990). 

2. P. Sikivie, Phys. Rev. Lett. 51, 1415 (1983). 

3. K. van Bibber et ah, Proposal for an Cosmic Axion Search Experiment at Lawrence 
Livermore National Lab (1990). 

4. J. Preskill, M. Wise, and F. Wilczek, Phys. Lett. B 120, 127 (1983); L. Abbott and 
P. Sikivie, Phys. Lett. B 120, 133 (1983); M. Dine and W. Fischler, Phys. Lett. B 
120, 137 (1983); M.S. Turner, Phys. Rev. D 33, 889 (1986). 

5. See, e.g., J. Caldwell and J.P. Ostriker, Astrophys. J. 251, 61 (1981); J.N. Bahcall 
and R.M. Soniera, Astrophys. J. Suppl. 55, 67 (1984); G.R. Knapp, S.D. Tremaine, 
and J.E. Gunn, Astron. J. 83, 1585 (1978); J.M.A. Danby and G.L. Camm, Mon. 
Not. r. Astron. Soc. 117, 50 (1957); J.M.A. Danby and T.A. Bray, Astrophys. J. 72, 
219 (1967); C.W. Allen, Astrophysical Quantities (The Athlone Press, London, 1973). 
For an excellent and clear discussion of the sun’s motion within the galaxy see D. Mi- 
halas and J. Binney, Galactic Astronomy: Structure and Kinematics (W.H. Freeman, 
San Francisco, 1981), Chapter 6. For practical assistance in transforming between 
the various astronomical coordinate systems — equatorial, ecliptic, and galactic — see 
P. Duffett-Smith, Practical Astronomy with your Calculator (Cambridge Univ. Press, 
Cambridge, 1981). 

6. K. Griest, Phys. Rev. D 37, 2703 (1988). 

7. For an elementary treatment of orbital mechanics, see, e.g., J.M.A. Danby, Fundamen- 
tal of Celestial Mechanics (Macmillan Co., New York, 1970), especially pp. 125-128, 
152-155, and 337-338. 


7 


8. A.J. Drukier, K. Freese, and D.N. Spergel, Phys. Rev. D 33, 3495 (1986); K. Freese, 
J.A. Frieman, and A. Gould, Phys. Rev. D 37, 3388 (1988). 

9. The importance of the microwave line due to axion-photon conversion as a probe of 

the galactic halo has also been emphasized by L. Krauss et al., Phys. Rev. Lett. 55, 
1797 (1985). ' 



FIGURE CAPTION 

Fig. 1. Local phase-space distribution function df /du for cosmic axions. The axion kinetic 
energy E a = u(m a v 2 / 2), which corresponds to a frequency of 980u (m a /10 -5 eV) Hz. 
The distribution function in the rest frame of the galaxy (tallest and narrowest), in 
the laboratory frame in June (shortest and broadest), and in the laboratory frame in 
December are shown. 
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